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Abstract 

We compute the third-order correction to the S-wave quarkonium wave functions 
|V>n(0)| 2 at the origin from non-Coulomb potentials in the effective non-relativistic 
Lagrangian. Together with previous results on the Coulomb correction and the 
ultrasoft correction computed in a companion paper, this completes the third- 
order calculation up to a few unknown matching coefficients. Numerical estimates 
of the new correction for bottomonium and toponium are given. 



1 Introduction 



The non-relativistic bound-state problem has a long history since the birth of quantum 
mechanics. Its systematic derivation from the relativistic quantum field theory of elec- 
trodynamics or chromodynamics (QCD) was developed more recently. Non-relativistic 
effective field theories [Tf2f3] together with dimensional regularization and diagrammatic 
expansion methods [I] now allow calculations of higher-order perturbative corrections, 
including all "relativistic" effects, to the leading-order bound-state properties, given by 
the solution of the Schrodinger equation. This is of interest in QCD for the lowest 
bottomonium state and top-antitop production near threshold, where non-perturbative 
long-distance effects can be argued to be sub-dominant, but perturbative corrections are 
large. 

The iS-wave energy levels are currently known at next-to-next-to-next-to-leading or- 
der (NNNLO)0 [5ll6f7f8] . except for the three-loop coefficient of the colour- Coulomb po- 
tential, but the corresponding wave functions at the origin, which are related to elec- 
tromagnetic decay and production of these states are completely known only at next- 
to-next-to-leading order (NNLO) [9|10|lllj . There exist partial results for logarithmic 
effects at NNNLO [T2|13lll4|ri5j . which can be related to certain anomalous dimensions 
and lower-order quantities. In [7j we computed the third-order corrections to S*-wave 
wave function at the origin from all terms in the heavy-quark potential related only to 
the Coulomb potential. In this paper we compute the contribution from the remaining 
potentials. A companion paper [16] deals with the Lamb-shift like contribution from 
ultrasoft gluons, thus completing the calculation of all bound-state effects at NNNLO, 
except for a few unknown matching coefficients. Our result is provided in such a form 
that these coefficients can be easily inserted, once they are computed. 

In contrast to the Coulomb corrections the calculation of the more singular non- 
Coulomb potential corrections leads to divergences, both in the calculation of the po- 
tentials themselves as in the insertions of these potentials in the calculation of the wave 
function at the origin. We employ dimensional regularization with d = 4 — 2e throughout, 
and provide a precise definition of all quantities, which corresponds to the MS subtrac- 
tion scheme. The technical details of this calculation together with an extension to the 
full S-w&ve Green function will be given elsewhere. 

2 Relating the leptonic quarkonium decay constant 
to the wave function at the origin 

We consider the two-point function 

(VV - g^q 2 ) I% 2 ) = i J d d xe^ (n\T(f(x)f(0))\ty, (1) 

1 Non-relativistic perturbation theory is an expansion in a s and the non-relativistic velocity v, while 
counting a s /v ~ 1, which implies a summation of the series in a s even at LO. We do not sum logarithms 
of a s lnu. 
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of the electromagnetic heavy-quark current = Q^^Q, choosing = (2m + E, 0) 
with m the pole mass of the heavy quark. The two-point function exhibits the S*-wave 
bound-state poles at E n , near which 

R( 2) E^_E n _Nc Zn , x 

Here N c = 3 denotes the number of colours. The residue Z n is related to the leptonic 
decay width r([QQ] n - ► l + l~) of the nth 5-wave quarkonium state by 

T([QQ} n - /+/") = y^T^i (3) 

with eQ the electric charge of the heavy quark in units of the positron charge, and a the 
fine-structure constant. Although there are no toponium states, and the cross section 
of top-antitop production is determined by the full two-point function, the residue Z n 
for n = 1 provides an approximation to the height of the broad resonance in this cross 
section. 

The electromagnetic current is expressed in terms of the non-relativistic heavy 
quark (z/>) and anti-quark (x) field operators via 

3 l = c v + D 2 X + . . . , (4) 

where the hard matching coefficients have expansions c„ = 1 + J2n ^ n ' ) ( a s/4vr) n , and the 
d v — 1 + d^(a s /4:Tr) + . . .. The central quantity in this paper is the two-point function 



G(E) = m ^ d — j d d xe iE *° (n\T( [tfo*x\(x) [xW](o)M 

|^(0) | 2 . . 

- En -E-te> [b) 

defined in non-relativistic QCD (NRQCD), whose poles define the wave functions at the 
origin and bound-state energy levels. At leading order, the wave functions and binding 
energies are given by |^°)(0)| 2 = (mC F a s ) 3 /(8irn 3 ) and E^ = —m(a s C F ) 2 /(An 2 ) 1 re- 
spectively (here and below C F = (N 2 - 1)/(2JV C ) = 4/3, C A = N c = 3). They receive 
perturbative corrections from higher-order heavy-quark potentials and dynamical gluon 
effects, hence E n = E® (l + Ek(a s /^) k e k ) and |^„.(0)| 2 = |^ 0) (0)| 2 (1 + £ fc (a s /4vr) fc / fc ). 
Using an equation-of-motion relation, we can replace D 2 in © by —mE, and we obtain 



c v 



E n ( \ , d v \ 

Cv 1 + — + ■ 

m V 6 1 



l^n(0)| 2 , (6) 



where terms beyond NNNLO are neglected. Inserting the perturbative expansions and 
defining Z n = \iJ n \0)\ 2 (1 + J2 k (u s /4ii) k z k ), results in 

z x = 2c« + f u (7) 
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4 16vr 2 M°) 



~2 



Z 3 



2c( 2 )+ C « 2 + 2 C « /l + /2 - 

3 maj. 

2cf +2c«( C i 2 ) + / 2 ) + (2 C ( 2 ) + C « 2 ) /l + / 3 



(8) 



167T 2 M°) 



ma 



(9) 



Note that e^, and depend on the principal quantum number n of the energy level, 
but we omitted a corresponding index to keep the notation short. The short-distance 
coefficients ci > , c\,> in the MS schema! are given in [T7lfT8] . The third-order coefficient 
<4 3 ) is not yet known except for the rif terms {rif is the number of quarks whose mass is 



set to zero) [19J. The one-loop correction to the derivative current can be obtained by 
applying a spin-triplet projection to the results given in [20], which gives 



-Cj, 



, u 16 
32 In — + — 

3 



m 



(10) 



Here the infrared- divergent part d^-' dlv ^ = — 16CV/e is subtracted and added back to 
the ultrasoft calculation, where it cancels an ultraviolet divergence [16]. The first- and 
second-order corrections e±, f\ and / 2 are summarized in [7]o The key quantity in the 
present work is the third-order correction f% = f% + f^ c + /3 s , which we split into three 
parts: f£ from the Coulomb potential given in [7] (see also [8]); f^ c accounting for 
all remaining potential insertions calculated below; and the ultrasoft correction /3 s = 

are finite, but all other expansion 



647r 2 5 us Vv calculated in [E]. Note that e x , fx, cW 
coefficients d!^\ f 2 and f%, f^ c , f£ s depend on a convention for regularizing 

(and subtracting) various ultraviolet and infrared divergences that arise in separating 
the contributions from the different scales. These divergences cancel in the expansion 
coefficients z^, since Z n is an observable. We regulate all ultraviolet (UV) as well as 
infrared (IR) divergences dimensionally, and adopt the MS subtraction scheme, but we 
also check explicitly that the pole parts of the various terms cancel in the sums ()9]) 
wherever possible. 



3 Potentials 

Starting from full QCD, potentials arise after integrating out hard and soft modes, 
and potential gluons [2T]. The potential NRQCD (PNRQCD) [2] effective Lagrangian 
contains a series of instantaneous interactions (potentials) with matching coefficients 
computed to an appropriate order in perturbation theory, and the interactions of ultrasoft 

2 The MS scheme is defined by the loop integration measure fM 2t d d k/ {2ir) d with jl 2 = /j, 2 e 7i3 /(47r) and 
subtraction of the pole parts in e. 

3 In Eq. (30) of the published version of [7], the term —En^/m in ijHJ) that comes from the expansion 
of q 2 around 4m 2 was incorrectly included in the formula for f£ . The term — 13/(8n 2 ) in Eq. (30) 
should therefore read — 15/(8n 2 ). This is corrected in [hep-ph/0501289^ 2]. 
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gluons. In this paper we do not consider the ultrasoft gluon interactions (see pSj), so 
the required Lagrangian is 



\ /.m Xm° / \ Am 



2m &m c 



+ / d 



^Vl (x + t) (V (r) + 5V(r, d) 



X 



x . 



(11) 



In the first line relativistic corrections to the kinetic energy are included; the second 
line shows the potentials with 5V the correction to the leading-order Coulomb potential 
Vq. The general form of the potential is colour- and spin-dependent. In writing (ITT)) we 
have made simplifications that apply to the computation of the Green function of colour- 
singlet, spin-triplet currents. We have also used that annihilation diagrams contribute 
to this Green function only from the fourth order. Since we work in dimensional regu- 
larization, the spin-triplet projection has to be carried out consistently in d dimensions, 
which implies tr(erV l ) = 2 {d — 1) for the Pauli matrices. Furthermore one must not use 
the e-symbol, and treat the commutator [cr%cr- ; ] as irreducible. For the calculation the 
momentum space representation of 5V is more convenient. It can be written in the form 
(q = p - p') 



5V 



Aita s Cp 
q 2 



V C - V 



l/m 



71-2 |q| + v 2 q 2 |V p : 

m m 2 v 2m 2 



P 



12 



(12) 



The strong coupling constant a s = a s (fi) is renormalized at the scale \i unless stated 
otherwise. The coefficients Vi can be expanded in a s , 



V; 



vf 0) + ^vf } + 



(2) 



0{al 



(13) 



with the superscript denoting the number of loops. The tree-level Coulomb potential 
V^- ^ is a leading effect and is therefore included in Vo(r) in the Lagrangian. For the 
calculation of the third-order non-Coulomb correction we need Vu m to two loops, and 
Vi/ m 2, V p to one loop. We also need the one-loop Coulomb potential to calculate the 

double insertion of Vq with Vy m2 and . 

The potentials must be determined in d dimensions, since the subsequent insertions 
are ultraviolet divergent. We now summarize the expansion coefficients. The tree- level 
coefficients are ITT] 



y(0) 



v (o) 

V l/m 



0. 



V (0) 



2e 2 



6-4e 



(14) 



The matching coefficients in PNRQCD receive contributions from hard momenta k ~ m 
and soft momenta k ~ ma s . Including counterterms (see below) the d- dimensional 
coefficients assume the form (/3 = HC^/3 — 4Tpri//3) 



v (i) 



7 + vq^ )ai(e) ' 



(15) 



4 



v (l) 

V l/m 



V (2) 
V l/m 



— ) (2C F C A + Ci 



+ 



2\2 



Wo 



2\2 



(16) 



(17) 



v (i) 



yd) 



7 11 

-O - — C A +Pov (e) 

3 6 



+ 



4-i 



2 y 



q 



I 1 



(18) 
(19) 



The calculation of the coefficients V^ m , V{^ 2 and V^ 1 -* results in IR divergences [6|22] . 
In writing the above equations, we subtracted the expression 



5V r 



c.t. 



a s C F 
6e 



C 3 A ^ + A(C 2 A + 2C A C F 



m\q\ 



16 C F 



, 16CU 

77T m z 



a. (P 2 + P /2 ) 



2q 2 



(20) 



which corresponds to a MS definition of the spin-projected potentials as PNRQCD 
matching coefficients. The counterterm ( 120]) is added back to the ultrasoft contribu- 
tion [16], where it cancels a corresponding UV divergence. The terms proportional to 
A)/e i n CE1) to (1191) are related to QCD charge renormalization. With the exception of 
61(e), which we need to second order in the e-expansion, all other coefficients are needed 
at the first order. Their expressions are 



a>i[e) 

h(e) 
b 2 (e) 



C A [11 -8e] —ATpUf 
9Z [i _ 2e] -C A [1- e 



e^r(l-6)r(2-e)r(e) fa 
(3-2e)T(2-2e) e' 

; ^r(|-6) 2 r(| + 6) 

7rir(l-2e) 



18 3 



C4C 



101 

^6~ 



4 

-ln2 

3 



r 2 

°A 



+ 



49 2 
36 Ca ~9 Cf . 



Tprij + eb 2 



(<:) 



^ 11 n (l,e) 

6^4 H Tpn f + e ?A ' ; 

3 27 27 * q 



(21) 
(22) 



(23) 
(24) 
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v^(e) = -^-fc A + ^T F + ev^\ (25) 

* )(e) = I^-i T ^ + e *' £) ' (26) 

where we checked the <i-dimensional expression for ai(e) in [23], 61(e) is from [TT], 6 2 (e = 
0) from |24j . and uf^(e = 0) for z = {g, m,p} are obtained from [6|25fi The 0(e) term 
of 62(e) is unknown; the corresponding terms v^ 1 '^ can be obtained from [25]. However, 
in the present work, we shall use this result only for a rough estimate. 

4 Third-order non-Coulomb correction to the wave 
function at the origin 

Since the relevant terms in the effective Lagrangian ffTTT) do not take us out of the 
quark- ant iquark sector of the Fock space, the calculation of G(E) can be mapped to a 
Hamiltonian problem in single-particle quantum mechanics after separating the trivial 
center-of-mass motion of the quark-antiquark pair. At leading order G(E) equals the 
Green function 

G (E) = (0|G |0> = (0| 1 . |0), (27) 

n.Q — £j — l€ 

of the Schrodinger operator Hq = V 2 /m — a s Cpjr with |0) denoting a relative position 
eigenstate with eigenvalue r = 0. The effect of the perturbation potentials 5V is taken 
into account by substituting H — > H = H + 5V, and then the Green function G(E) 
is systematically expanded in powers of a s . The third-order non-Coulomb correction 
corresponds to the expression 

5 3 G = 2 (0\G 5V 1 G 5V 2 G \0) - (0\G SV 3 G \0) (28) 

with 8V\ the one- loop Coulomb potential contribution to 5V, and 8V2, SV 3 the second 
and third order non-Coulomb potentials in momentum space given by 

5% = ^ (27F)d_1 6id ~ 1){P ~ p,) + ^ + ^ + ^ (0) ' (29) 

W* = WiJl + W$* + Wf\ (3°) 

where the first term in 8V2 is due to the relativistic kinetic energy correction in (fTTj) . 
and refers to the corresponding terms in (fl~5j) to (fT9|) including the prefactors in 

(|T2|) . The result of the potential insertions (|28|) is matched to the pole structure of 
G(E), which allows us to extract the bound-state energy and its wave function at the 



4 The coefficient u^ 1 ^ corresponds to —d^ in [6]. 
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origin. The technical details of the calculation will be discussed elsewhere. Here we 
briefly summarize the main results. 

The expression (]28|) should be considered as an infinite sum of loop diagrams defined 
in dimensional regularization. Since the (d — l)-dimensional Coulomb Green function 
Go is not known in closed form, one must separate the divergent subgraphs. Since these 
contain only a finite number of loops and since the divergent part is local, one calculates 
these subgraphs in d dimensions and factorizes the pole part from the remainder of 
the diagram. All non-divergent diagrams can be computed in d = 4 using the known 
three-dimensional Coulomb Green function. For the divergent part of (128]) we find 



6*G 



dm 



1 



7 

72 



Cl 



-C 2 



% CaCf + Po ( 2l 



36 




CaCi 




CfTf 



49 
216 



f 



\Ca + \c F 



a 2 C F 



(0|G ^G |0) 



IT 



(0|G |0> 



(31) 



with L m = ln(/i/m). In the residue Z n the pole part in the last line cancels against 
a corresponding term in the contribution 2c^fi in <$f), while the divergent part pro- 
portional to Gq(E) combines with similar terms in the ultrasoft contribution, and must 
cancel with 2c^ . From the finite part of S^G, we extract the non-Coulomb contribution 
to fz, and obtain 



fnC 
J3 

64vr 2 



qC f + —C A C F + -C 2 A C F + Po ( -C 2 F + 2C A C F \ 



cl 



3 14 7Si 
2 + M~~ 



A 

109 



CaC f 



226 
~27~ 



3n 



Q C A C F T F n f + f3o<C F 



16 10 



81n2 37 5 _37S 1 
3 3n 3n 2 6 



75 



3n 16n 2 



■K 2 n 45i 2nS 2 



9 



+ 



+ C a Cf 



15 



5 
n 



■K 2 n 



1 1 

-Cf + -C A C 2 



L m + 



1 

15 



^C 2 T F 



Lm + 



rXlC 

c tp,3 

64vr 2! 



2Si + nS 2 



12 3n 



L + 
Si 



-cl 



+ CaC f 



-CaC 2 f 



9 n 



Si 
2 



(32) 
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where L = ln(n/^/(mCVa s )), and the unwritten argument n of the harmonic sum Sn = 
S a (n) = Z)fc=il/^ a * s the principal quantum number n. The non- logarithmic partly is 
given by 



„nC 
C ip,3 

64vr 2 



+ 



137 49tt' 



36 



432 



25 
6n 



7061 50tt 2 



486 



SI 



+ 



V 27 
3407 5tt 2 



1475 
108n 

81n2 



7T 

9rz 



35 

12^2 + 

321 



3 14 
Sl 1 2 " 3^ 



In 2 



/353 



133 



h 37 + 1 
3n n 2 

187 



32n 2 V 54 

37 Si 



7S 1 
IT 

16 

3n 



7^2 
~ IT 

16 In 2 



C% 



9 



12 



37 
12 



So[ — + — 



2 

3n 



CaCf 



432 



9n 

41n2 



In 2 



8 81n2 



108 3n 



9 



AS 2 



+ 



361 
108 

3391 



18 ' 3 
49 In 2 109 



16 

3n 



+ 



108 

5tt 2 



486 648 
f 1027 19 



18n 

21n2 
27 

25 



4Si 
IT 

10951 
36 

118 



1 4 
15 + 15n 



25i 
l5~ 



ClT, 



F-LF 



C aC fTf Tlf 



125 59Si 



\ 648 



-St 
24 ' 



6n ^ 24n 2 



27n ' 24n 2 
35tt 2 



27 



108 



+ 



10 

J + 

917T 2 

144 



1 45 
3n IQn 2 



ii 2 n 



ll7T 2 n 57T 2 

27 + 16n + 
2n,S 2 



CpTpUf 
4nS 3 



2nS 2 ,i 



i 

An 



5ir 2 n 
24 



5-1 



1 

2n 



7r 2 n 



/3 5n 

S 2 [ — I 

\2 4 



nS^i + 2nS , 3 > CaC f 



where S a ,b = S at b(n) = Sfc=i Sb{k)/k a is a nested harmonic sum, and 



,(1,0 



+ 



,(1,0 
J? 

12 



_p 

12 



6 



C 2 - ^ 



2 , 



(33) 



(34) 



denotes the contribution from the unknown 0(e) terms of the potentials. Note that S e 
is independent of the principal quantum number n. 



5 For some of the terms in our calculation, we obtained analytic expressions which can be evaluated 
only for any given value of n. To obtain the presented result for general n, we made a rather general 
ansatz for the possible dependence on n such as harmonic sums and determined the rational coefficients 
of the various structures from the analytic result for n = 1, 2, . . . , 8. The correctness of this result was 
then checked for all further n up to n = 30. 
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5 Numerical results 



We briefly discuss the numerical significance of the non-Coulomb correction. The num- 
bers below should be considered as indicating the typical size of the correction. The 
correction itself is strongly scale-dependent (as well as scheme-dependent due to (1201)) 
and a physically relevant result is only obtained by combining all terms in (Q. This 
cannot be done, since and S e are still unknown. 

The numerical expressions of the non-Coulomb wave function correction to the lowest 
three states n = 1,2,3 are given by (Wn(0)| 2 c = \^\0)\ 2 (a s /(4ir)) 3 f 3 nC ) 



«i(0)lL 



Ss\M0)\la 



Ss\M0)\la 



( ma f F ^ A \ (149.3 - 6.9n/) L 2 + 0.9L 2 n + 3.5L L m 

+ (449.8 - 21.9n/) L + 0.8L m + (-149.7 - 3.1n/) + 8 € 
{ma C F f a| r g _ ^ + Q Q ^ + g ^ ^ 

64:71 7T L 

+ (211.7 - 13.5n/) L - AAL m + (-217.9 + 2.7n/) + 5 e 
(ma s Ci?) 3 a 3 



2167T 7T 

(89.7-8.8n/)L-6.7L 



149.3 - 6.9n/) L 2 + 0.9L 2 m + 3.5L L 
(-218.3 + 5.1n/) + 6 e 



(35) 



(36) 



(37) 



Since some of the 0(e) parts of the potentials, 5 t , are not known, we perform a crude 
estimate. We derive the expressions for \ i = m,q,p from [25], and estimate 
&2 /^(e = 0) by the corresponding ratio for &i(e), allowing for a factor ±2 variation 
in this ratio. In this way we find that 5 e lies between —3 and +8, is dominated by 
but is more than an order of magnitude smaller than the other non-logarithmic terms. 
Below, we therefore set 5 e to zero. The potential scale is of order ma s Cp, while the 
natural choice for the hard factorization scale is m. In the following we show the size of 
the non-Coulomb correction for two scale choices, \xb = mC F a s (fiB) and /i = m. 

"Toponium" wave function at the origin. For the top-antitop production cross section 
near threshold, the height of the smeared resonance is related to the wave function at 
the origin of the n = 1 ground state. The two scales discussed above correspond to 



cvi n/ " 5) (M = 0.140 (ji B 
(135|) . we obtain 



32.6 GeV) and a 



(«/=5), 



m t 



0.107 (m t = 175 GeV). From 



«i(0)|£c a20* 



7T 



l4 0) (o)| 2 



■165.1 + 0.8 In (a s C F ) + 0.9 \n 2 (a s C F ) 



-0.14 



l^ 0) (o)| 2 



7T 



-165.1 + 340.3 In 



o> s Cf 



+ 114.7 In" 



a s C F 



(38) 
0.36, (39) 
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where \i — (a* — m t) is used in the first (second) line. Our estimate for 5 t affects 
these numbers by less than 0.01. For comparison, the NNNLO Coulomb correction is 
( ^3|V ; i(O)|c/l'0i O ' ) (O)| 2 = —0.04 with \i — v — \i B (see [7J). The result above indicates a 
potentially very large third-order contribution from the non-Coulomb potentials. A sim- 
ilar observation is made for the ultrasoft contribution [IB] , and since the non- logarithmic 
terms of the non-Coulomb and ultrasoft contribution have opposite sign, a cancellation 
is possible. 



Bottomonium. Here the QCD corrections to the wave function at the origin are known 
to be rather large, and the application of the perturbative calculation is certainly ques- 
tionable for excited states [7J. Therefore, mainly the bottomonium ground state T(IS') is 
of phenomenological interest. The bound-state scale corresponds to oT^^^b) = 0.30 
(fi B = 2GeV), and the hard scale to ai™ /=4) (5 GeV) = 0.21 (m b = 5GeV). The non- 
Coulomb correction to the Y(IS') wave function at the origin is 

Willie _ * 3 sM ( _ 162 + g HasCp) + i^ {as c F ) ) = -1.4, (40) 



l^ 0) (o)l 2 



, -162.0 + 362.2 In — - + 121.6 In 2 — - =1.5. (41) 
|^i 0) (0)| 2 TT V \a,C F J \a s C F JJ 

Once again, the contribution of 5 e does not seem to be important. For comparison, 
the NNNLO Coulomb correction is ^^(O^/l^f^O)! 2 = -0.47 for n = fi B - In view 
of the size of this third-order correction the validity of a perturbative treatment of the 
leptonic T(IS') width relies on cancellations that might occur in the complete third-order 
correction. 



6 Summary 

We computed the third-order correction to the S'-wave quarkonium wave functions at the 
origin originating from insertions involving potentials other than the Coulomb potential. 
The UV as well as IR divergences are regulated dimensionally in a way that is consis- 
tent with other parts of the calculation, performed in the context of the diagrammatic 
threshold expansion. Together with the previous calculation of the Coulomb potential 
contributions [7115] and the ultrasoft contribution [TB] this completes the bound-state part 
of the third-order calculation. The remaining unknowns are related to the finite parts 
of three-loop matching coefficients (03, c^) and the O(e) parts of one- and two- loop 
coefficients (parameterized here by 5 e ). Numerical estimates of the new non-Coulomb 
contribution suggest important third-order effects even for top quarks, but a definitive 
conclusion can only be drawn once the missing matching coefficients are computed. 
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